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Abstract. The (if, T)-phase diagram of the multiferroic perovskite TbMn03 was 
studied by high- resolution thermal expansion a(T) and magnetostriction AL(H)/L 
measurements. Below T/v — 42 K, TbMn03 shows antiferromagnetic order, which 
changes at Tfe — 28 K where simultaneously a spontaneous polarization P\\c develops. 
Sufficiently high magnetic fields applied along a or 6 induce a polarization flop to P\\a. 
We find that all of these transitions are strongly coupled to the lattice parameters. 
Thus, our data allow for a precise determination of the phase boundaries and also yield 
information about their uniaxial pressure dependencies. The strongly hysteretic phase 
boundary to the ferroelectric phase with P\ \a is derived in detail. Contrary to previous 
reports, we find that even in high magnetic fields there are no direct transitions from 
this phase to the paraelectric phase. We also determine the various phase boundaries 
in the low-temperature region related to complex reordering transitions of the Tb 
moments. 
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1. Introduction 

The coupling between magnetic and ferroelectric order parameters in so-called 
magnetoelectric multiferroic materials is of great current interest [1-5]. After the 
rediscovery of this mechanism a few years ago, an intense search has started for 
materials with strong coupling of the spontaneous polarization P and the spontaneous 
magnetization M [6-8]. This search is enhanced, on the one hand, by the demand for 
new promising components for device design [9]. On the other hand, the fundamental 
aspects of the magnetoelectric coupling in many systems are far from being understood. 
Thus, a detailed determination of the complex and rich phase diagrams observed 
in multiferroic materials is important. Considering the multiferroic orthorhombic 
manganites i?MnO"3 (R = Gd, Tb, Dy), most of the previous investigations concentrated 
on the ordering phenomena related to the Mn ions, while less is known about the 
ordering of the rare earth ions in magnetic fields [10-12]. Recently, we have shown 
that measurements of thermal expansion and magnetostriction by high-resolution 
dilatometry are a powerful method to investigate the temperature and magnetic-field 
phase diagrams of multiferroics, because both the magnetic and the ferroelectric order 
strongly couple to the lattice parameters [13]. These couplings reflect pronounced 
uniaxial pressure dependencies of the respective transition fields and temperatures, 
which can be analyzed by Clausius-Clapeyron and Ehrenfest relations for first- and 
second-order phase transitions, respectively [14]. 

Below the Neel temperature T N ~ 42 K the Mn spins of TbMn0 3 
develop an incommensurate sinusoidal antiferromagnetic alignment with wave vector 
(0,kMn(T),0) [15] (using the Pbnm setting of the orthorhombic unit cell). According 
to Harris et al. [16], this phase is called the high-temperature incommensurate 
antiferromagnetic phase (HTI). At T FE ~ 28 K another transition occurs, leading also 
to an incommensurate, but cycloidal ordering of the Mn moments along (0, kMn(T), 0). 
Because the latter transition breaks the inversion symmetry of the crystal, a spontaneous 
polarization P\ \c can appear below Tfe even in zero field [17,18]. This phase is called the 
low-temperature incommensurate AFM phase (LTI). In zero field, an incommensurate 
AFM ordering of the Tb ions has been proposed [11], which sets in below T^ h = 7 K. 
The above-described ordering phenomena of the Mn ions hardly change for magnetic 
fields H\\i < 4.5 T (i — a, b, c). However, in higher magnetic fields significant changes 
in the ferroelectric LTI phase occur: sufficiently high fields (H > i?FE,c) parallel to 
the a or b axis induce a polarization flop from P\\c to P\\a. The critical field strengths 
Hfe,c depend on their orientation and on temperature. In general, i?FE,c is smaller 
for H\\b (# FEC > 4.5 T) than for H\\a (#fe,c > 9.5 T) [12]. The polarization flop is 
accompanied by a change in the modulation wave vector, which becomes commensurate 
for H > i?FE,c (LTC phase) [15]. If a magnetic field H > 8 T is applied parallel to the 
c axis, the sample is forced into a canted antiferromagnetic ordering (cAFM phase) and 
the ferroelectric order is suppressed completely [12]. 

In this paper we present a study of the linear thermal expansion coefficient 
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Oii(T) = j- 9L q^ and the linear magnetostriction AL ^ H ^ = [Li{H) — Lj(0)]/Lj(0) 
of TbMnOs. Here, denotes the length of the sample parallel to the different 
crystallographic directions i — a, b, and c. The TbMnOs single crystal was cut from 
a larger crystal grown by floating-zone melting in an image furnace [12]. The paper 
is organized as follows: in the next section we present zero field measurements of the 
linear thermal expansion coefficients. Section 3 deals with the phase transitions which 
are attributed to the Mn sublattice. We performed measurements along all three lattice 
directions in longitudinal magnetic fields up to H = 16 T, i.e., in all cases the field 
was applied parallel to the measured crystal axis (H\\i with i — a, b, and c). The 
phase diagram derived from our data and the uniaxial pressure dependencies of various 
phase boundaries are discussed in the fourth section. Section 5 considers the ordering 
transitions related to the Tb moments and the last section gives a summary. 



2. Linear thermal expansion in zero field 

The thermal expansion coefficients oti(T) of TbMnOs measured in zero magnetic field are 
shown in figure [U As may be expected from the orthorhombic structure, the TbMnOs 
single crystal shows a strongly anisotropic thermal expansion. Several anomalies are 
detected. These anomalies can be attributed to the different phase transitions according 
to previous publications [15]. The sharp anomalies at 41.5 K signal the Neel transition 
of the Mn ions. At Tfe = 27.6 K TbMnOs undergoes a second phase transition from the 
HTI to the LTI phase, accompanied by a pronounced anomaly in at(T) for i = a and 6, 
whereas the effect in a c (T) is much smaller. The shape of the anomalies at Tn and T F e 
is typical for second-order phase transitions. The uniaxial pressure dependencies of T c 
are described by the Ehrenfest relation 
dT c _ Aa t 

-d^ ~ C ~Ac ■ (1) 
Here, V m is the molar volume, A«j (Ac) denotes a jump in the thermal expansion 
coefficient (specific heat) and i is the measurement direction of and the direction of 
uniaxial pressure pi. Because any ordering transition causes a decrease of the entropy, 
the specific heat anomaly is always positive (Ac > 0) and the sign of dT c /dpi is given 
by the sign of Actj. Thus, T N shifts to higher temperature for uniaxial pressure Pi along 
the b or c axis and decreases for p a . The changes of T F e for p a and pb have the same signs 
as those of T N , i.e., T FE decreases for p a and increases for p&. Pressure parallel to the c 
axis has only a minor effect on Tfe, since practically no change Aa c is observed at T F e- 
The ordering of the Tb sublattice at T^ h = 7 K also causes pronounced anomalies for 
all three crystallographic directions. From the signs of the thermal expansion anomalies 
we conclude that the Tb ordering is stabilized for uniaxial pressure along a or 6, while 
it is suppressed for p c . 

Apart from the anomalies related to the different transition temperatures there are 
additional anomalous features visible for all three directions. Firstly, kink-like anomalies 
appear at T kink ~ 34 K. We suspect that these kinks may be related to the slope change 
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Figure 1: Thermal expansion of TbMnC>3 measured along the a, b, and c axis. 



in the temperature dependence of kMn(T) reported by Kajimoto et al, although the 
latter has been observed at a somewhat larger temperature T ~ 38 K [11]. Secondly, 
a a shows a pronounced broad minimum (maximum) centered around 80 K and a 
weaker minimum is also present in a c . Such broad extrema are typical indications for 
Schottky contributions arising from a thermal population of low-lying excited states. 
The fact that this thermal population contributes to the thermal expansion coefficients 
reflects strong uniaxial pressure dependencies of the relevant energy splitting(s) between 
the ground state and the excited state(s). Again, the signs of the uniaxial pressure 
dependencies are given by the signs of the anomalous «j [19-21]. In TbMnC>3 these 
Schottky contributions most probably arise from the 4/ multiplet of the Tb 3+ ions. 
The Hund's rule ground state of the free Tb 3+ ion has a total orbital momentum J = 6 
and is 13-fold degenerate. In a crystal field of orthorhombic symmetry this degeneracy 
is completely lifted and it is obvious that the energy splittings between the 13 singlet 
states may strongly change with pressure, since the crystal field will depend on pressure. 
A more quantitative analysis of the Schottky contributions to is not possible at the 
present stage, since it would require a detailed knowledge of the different singlet states 
and of their energy splitting. 

3. Transitions related to the Mn ordering 

We will divide the discussion of our data in two parts. In this section we will concentrate 
on the phase transitions related to the Mn sublattice. The ordering of the Tb ions is 
discussed in section [5] This division does not mean that the Mn and the Tb sublattice 
of TbMnC>3 act independently from each other, but it is reasonable to assume that the 
coupling between the two sublattices is not too strong. 
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Figure 2: Panel (a): Field-cooled (FC) thermal expansion a a {T) for H\\a measured 
with increasing temperature. For clarity the curves for different fields are shifted by 
—5 • 10~ 6 /K with respect to each other. We find a pronounced hysteresis around T FE C . 
This is exemplified in the inset by the a a (T) curves measured in a field of 10 T with 
increasing and decreasing temperature. This hysteresis is also seen in magnetostriction 
measurements Aa(H)/a recorded with increasing (solid lines) and decreasing (dashed 
lines) field shown in panel (b). Panel (c) displays an enlarged view on the high-field 
a a (T) anomalies, which both have additional shoulders on their high-temperature sides, 
indicating the presence of two transition temperatures Tfe,c and Tfe- 



3.1. Measurements in magnetic fields H\\a 

Measurements of a a (T) in a longitudinal magnetic field H\\a up to 14 T are shown in 
figure [2] (a). Above about 10 K the a a (T) curves in 4 T and in zero field are quite 
similar. For higher magnetic fields H > H EEC ~ 9.5 T an additional anomaly appears 
at Tfe,c which signals the phase transition from the LTI (P||c) to the LTC phase with 
P\\a. This first-order transition shows a broad hysteresis. The hysteretic behavior is 
presented in the inset of figure [2] (a) by measurements of a a (T) in 10 T with increasing 
(dT/dt > 0) and decreasing (dT/dt < 0) temperature. The LTI-to-LTC transition 
can also be detected by measurements of the magnetostriction, i.e. the magnetic-field 
induced length change at a constant temperature. In figure [2] (b) magnetostriction 
measurements at T = 5 K and 20 K are presented. In both curves a jump-like expansion 
of Aa(H)/a is observed as the phase boundary between the LTI and LTC phase is crossed 
as a function of increasing H. For decreasing H, the anomalous expansion is reversed 
and a jump-like contraction signals the LTC-to-LTI transition. As in a a (T), this first- 
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order transition shows a hysteresis, which is strongly enhanced at lower temperature. 

As is also seen in figure [2] (a), a a (T) in if = 14 T shows a pronounced anomaly 
at ~ 28.5 K. At first glance, the peak-like shape and a small hysteresis (not shown) 
indicate that there is one first-order phase transition at this temperature. However, a 
closer inspection of the high-field anomalies for if = 12 T and 14 T reveals that both 
anomalies are asymmetric with additional shoulders on their high-temperature sides; 
see figure [2] (c). This observation indicates that in high magnetic fields two separate 
transitions have to be distinguished. A natural interpretation is that with increasing 
temperature a first-order transition from the LTC to the LTI phase takes place and 
this transition is followed by a second-order transition from the LTI to the HTI phase 
at a slightly higher temperature (~ 0.5 — 1 K). In this respect, there is no qualitative 
difference between the sequence of transitions in intermediate (< 10 T) and in higher 
fields (> 11 T). This conclusion is in contrast to previous publications [10,15] which 
proposed direct transitions from the HTI to the LTC phase in TbMnOs for H > 11 T. 

The Neel temperature Tn shows only a weak increase of about 2 K in the field range 
up to 14 T, i.e., the HTI phase is slightly stabilized with increasing H\\a. The field 
dependence of T kink is closely linked to the critical field ifp E c ~ 9.5 T. For H < if EE c , 
T kink ~ 33.5 K is nearly field independent and jumps to T kink ~ 28.5 K at if EE C J3 

3.2. Measurements in magnetic fields H\\b 

The thermal expansion coefficient a&(T) is presented in figure [3] (a) for some 
representative magnetic fields H\\b. Similar to the measurements in H\\a, the Neel 
temperature weakly increases with field and the anomaly at T kink shows a jump-like 
decrease at the critical field if FE c ~ 5 T. In a magnetic field of 16 T, a pronounced 
anomaly at T FE C = 26.5 K signals the first-order LTI-to-LTC transition. Again this 
anomaly has an additional shoulder on the high-temperature side, which arises from the 
second-order LTI-to-HTI transition at T EE ~ 28 K. For if = 8 T these two transitions 
are well separated from each other, T EE ~ 26 K and T E e,c — 20 K, and their different 
order is reflected in different shapes of the respective anomalies. In figure [3] (b) we 
compare the anomalies of cn&(T) at T FEj c for H = 8 T and H = 16 T by shifting the 8 T 
curve on top of the 16 T curve. No further scaling is applied. The additional shoulder 
of the 16 T curve is clearly seen and a subtraction of the shifted 8 T curve even allows 
for a quantitative estimate of its magnitude. As shown in the inset of figure [3] (b), the 
additional shoulder in 16 T is of comparable magnitude as the anomalies due to the 
HTI-to-LTI transitions at T EE observed in lower fields. Thus, our data for H\\b allow for 
the same conclusion as drawn above from our data for H\\a. For both field directions 
we do not observe direct HTI-to-LTC transitions. 

Figure 0] (a) displays the anomalies of oii(T) at T FE c measured with increasing 
and decreasing temperature in H = 10 T. The shape of the anomalies is typical for 

| The kink-like anomaly is only present in measurements of oti(T) (i = a, b, c) with increasing 
temperature. If ctj is measured with decreasing temperature, no kink appears. 
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Figure 3: Panel (a): FC thermal expansion a:&(T) for H\\b measured with increasing 
temperature. For clarity the curves for different fields are shifted by 4 ■ 10~ 6 /K with 
respect to each other. Panel (b) compares the anomalies of ctb{T) at T F e,c measured 
in H = 8 T (left and bottom scales) and in H = 16 T (right and top scales). The 
additional shoulder on the high-temperature side of the anomaly in 16 T is obvious. As 
is shown in the inset, the magnitude of this shoulder is comparable to the anomalies 
due to the HTI-to-LTI transitions at Tfe observed in lower fields. 



a first-order transition and there is a clear hysteresis. The different peak heights are 
not too surprising because this simply means that not only the value of T F e,c but also 
the width of the transition depends on the sign of the temperature drift. A much 
more surprising feature is found for lower fields, as is shown for H = 6 T in panel (b). 
Here, we observe a double-peak structure when the phase boundary from the LTC 
to the LTI phase is crossed with increasing temperature. However, for the opposite 
direction, i.e. dT/dt < 0, a single peak signals the LTI-to-LTC transition. The same 
anomalous behavior is observed in our magnetostriction measurements. As shown in 
figured] (c), a jump-like contraction of the b axis at H ~ 5.4 T signals the LTI-to-LTC 
transition as a function of increasing H at T = 10 K. This is typical for a first-order 
transition and the corresponding field derivative of Ab(H)/b shows a peak of negative 
sign; see figured] (d). On decreasing the field again, we do not, however, observe the 
expected jump-like expansion at the LTC-to-LTI transition. Instead there is even a 
further decrease of Ab(H)/b in a restricted field range, which causes a double-peak in 
the corresponding field derivative. These highly anomalous double-peak structures only 
appear at the LTC-to-LTI transitions at temperatures below about 20 K and for H\\b. 
In other compounds, double peaks of a(T) have been observed on heating through 
glass-like transitions and it has been found that these peaks depend on the rate of the 
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Figure 4: Hysteresis of the transition between the LTI and the LTC phases for H\\b. 
Panel (a) and (b) show a&(T) measured with increasing (•-) and decreasing (o • • •) T in 
a field of 10 T and of 6 T, respectively. In 6 T, the anomalies are of completely different 
shapes for different signs of dT/dt, while only the widths of the anomalies change in 
H — 10 T. Panel (c) shows magnetostriction data recorded with increasing ( — ) and 
decreasing (• • •) magnetic field at different T. For different signs of the field sweep 
the anomalies only weakly change at 25 K, while they completely change their shapes 
for T < 20 K. Panel (d) displays the field derivatives of the magnetostriction data at 
T = 10 K for increasing (•-) and decreasing (o • • •) field. The LTI-to-LTC transition 
causes a single peak, while a double peak is observed at the LTC-to-LTI transition, in 
complete analogy to the anomalies of a&(T) shown in panel (b). 

previous cooling process [23,24]. In order to test whether such memory effects also exist 
in TbMnOs, we varied the rate of the increasing field between 20 mT/min and 2 T/min, 
but we did not observe any influence in the subsequent measurements with decreasing 
field. At present, the origin of the anomalous double-peak structures remains unclear. 

3.3. Measurements in magnetic fields H\\c 

According to Ref. [15], no additional anomalies arising from the Mn sublattice are 
expected in magnetic fields along the c axis up to H ~ 7 T. Larger fields induce a 




Figure 5: Panel (a): FC thermal expansion a c (T) for H\\c measured with increasing 
temperature. For clarity the curves for different fields are shifted by 4 ■ 10~ 6 /K with 
respect to each other. Panel (b) compares the a c (T) curves (not shifted) in the field 
range between 4 and 7 T and the inset shows an expanded view of the anomalies at T N . 



canted AFM ordering of the Mn spins and the ferroelectric polarization is suppressed. 
Figure El (a) shows the thermal expansion coefficient a c (T) in fields up to 4 T and, 
indeed, these curves do hardly change above about 12 K. The drastic changes in the 
lower temperature region arise from the Tb sublattice and will be discussed in section El 
Figure [5] (b) presents the a c {T) measurements in the field range from 4 T and 7 T. In 
6 T, a c (T) already shows a broad minimum around 15 K, which changes into a well- 
pronounced anomaly at ~ 16.5 K in the 7 T curve. We attribute this anomaly to the 
transition from the LTI to the cAFM phase. Unfortunately, the investigation of this 
phase boundary by high-resolution dilatometry was not possible at higher fields for H\\c 
because of strong torque effects. Usually, the sample is clamped in the dilatometer by 
only a very small pressure, which was, however, not sufficient to prevent a rotation of the 
sample when the phase boundary of the cAFM phase was crossed. Thus, a fitting was 
constructed to safely fix the sample orientation, but the torque effects were sufficiently 
strong to crack the sample in this setup. 

4. (H, T)-phase diagram and pressure dependencies 

From our thermal expansion and magnetostriction measurements we obtain the phase 
diagram shown in figure El The phase boundary between the LTI and the LTC phase, 
measured as a function of increasing and decreasing field H\\a and temperature, has 
been determined for the first time. We find a very strong hysteretic behavior, especially 
at low temperatures. Concerning the phase boundaries between the paramagnetic and 
the HTI phases our data well agree to previous investigations [10, 12, 15] for all three 
field directions. In the field range below about 10 Tesla, this is also the case for the 
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Figure 6: Phase diagram of TbMnOs based on the measurements of aj(T) (triangles) 
and ALi(H) / Li (squares) measured along i = a, 6, c with H\\i. Filled and open symbols 
are obtained with increasing and decreasing temperature and field, respectively. The 
dashed lines depict the high-field phase boundaries obtained in Ref. [15], because we 
could measure our sample only up to 7 T for H\\c (see text). 



New features in the phase diagram of TbMnOz 



11 



boundaries between the HTI and LTI phases. In contrast to the earlier publications, 
however, our data clearly show that there are no direct transitions from the HTI to the 
LTC phase in higher fields, neither for H\\a nor for H\\b. Instead of a direct transition 
from the HTI to the LTC phase as a function of decreasing temperature, the system 
seems to pass always through the LTI phase, before the LTC phase is entered. This 
observation is of importance for microscopic models describing the symmetry changes 
at the various phase transitions. 

Besides the bare position of the different phase boundaries, the thermal expansion 
data also yield information about their uniaxial pressure dependencies (see eq. (pQ)). 
As discussed in detail in Ref. [14], the opposite signs of the anomalies of a a and ab 
at both, Tn and Tfe, mean that both transitions depend on the degree of distortion 
of the GdFe03-type structure, which can be parameterized by the magnitude of the 
orthorhombic splitting e = (b — a)/(b + a). The main idea is that the ferromagnetic 
nearest-neighbor exchange J™ in the ab planes is weakened with increasing e, as it 
is the case in the i?Mn03 series for a decreasing ionic radius from R — La . . . Dy. In 
addition, the anisotropy of the antiferromagnetic next-nearest-neighbor coupling J^nn 
increases, since J^™ increases along b and decreases along a. As a consequence, the 
zero field magnetic ground state changes from an A-type AFM for R = La . . . Gd to 
an E-type AFM for the HoMn03 in perovskite structure^ Depending on temperature, 
magnetic field and the radius of R, the competition of J™ and J^nn can a l so cause 
incommensurate AFM structures for R = Eu...Ho [27,28]. Since uniaxial pressure 
along a (b) will increase (decrease) e, the corresponding uniaxial pressure dependencies 
of T N of TbMnOs can be straightforwardly traced back to a decreasing (increasing) 
J™ due to the changes of the orthorhombic splitting. This result is identical to our 
conclusions concerning the uniaxial pressure dependencies of T N of GdMnOs for p a and 
Pb. For uniaxial pressure p c , we find an increase of T N in both compounds. As discussed 
in Ref. [14], the GdFe03-type distortion is characterized not only by the magnitude of s, 
but also by a decreasing lattice parameter c. Thus, one might expect a negative dT N /dp, 
arising from a partial increase of the GdFe03-type distortion lowering J™- However, 
T N does not only depend on the couplings J™ and J^nn acting within the ab planes. 
A three-dimensional ordering requires a finite coupling J^ FM along the c direction, and 
the positive <9Tn / dp c suggests an increase of J^ FM under pressure along c. 

Concerning the HTI-to-LTI phase boundary and its pressure dependencies, the 
phenomenology of TbMnOs is not directly comparable to our previous results on 
GdMnOs [13,14]. In TbMnOs, the transition to the LTI phase is accompanied by a 
finite polarization already in zero magnetic field, whereas in GdMnOs a finite field H\\b 
is necessary to induce ferroelectricity. In addition, the ferroelectric phase of GdMnOs is 
not entered directly from the the paraelectric high-temperature incommensurate phase 
(termed ICAFM in Refs. [13, 14]) but from a canted AFM phasejjj] Concerning the 

§ -RM11O3 with rare earth ions smaller than Dy usually crystallize in a hexagonal structure [25,26] 
|| Note that the magnetic structures of GdMn03 have not yet been determined unambiguously. The 
proposed phases are based on the observed weak ferromagnetism [10, 29] and on X-ray diffraction 
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uniaxial pressure dependencies for p a and Pb, we observed a clear anti-correlation in 
GdMnOs between Tn and T c (the latter signals the HTI-to-cAFM boundary) on the one 
hand side and T FE on the other. This anti-correlation suggests a competition between 
colinear spin structures, either incommensurate or A-type AFM, and non-colinear spin 
structures, which allow for additional ferroelectricity [17,18]. In contrast to GdMnOs, 
the uniaxial pressure dependencies of T N and T FE of TbMnOs have the same signs. An 
increase of the orthorhombic splitting by pressure would decrease both, Tn and TpE- 
This is in agreement with the observed lower values of both, Tn and Tfe of DyMnOs 
where the orthorhombic splitting is larger than in TbMnOs [27,28]. Our finding suggests 
that the orthorhombic splitting of TbMnOs is already larger than the optimum value of e 
needed to establish a ferroelectric ordering with a maximum value of T F e- Summarizing 
the above discussion, we conclude that the competition between a multiferroic and an 
A-type AFM ground state dominates in GdMnOs, while in TbMnOs the dominant 
competition is between the multiferroic and an E-type AFM ground state. 

Apart from the phase boundaries related to the Mn subsystem, figure[6]also contains 
various phase boundaries in the low-temperature range arising from transitions of the 
Tb ions. The corresponding measurements, which have been used to trace these phase 
boundaries will be discussed in the following section. We find evidence that only 
one ordered phase of the Tb moments exists for H\\a, while for H\\b it is possible 
to distinguish at least two different phases. For H\\c the situation is more complex 
because a clear separation of different phases is not possible from our data. 

5. Transitions related to the Tb ions 

Even though the Tb ions have a strong magnetic moment which may interact with 
the Mn moments, no systematic investigation on the field dependence of the low- 
temperature phase transition has been published so far. According to Refs. [11,12] an 
incommensurate AFM ordering of the Tb moments occurs in zero field below T^ h = 7 K. 
In Ref. [10] Kimura et al. present several phase boundaries in the relevant temperature 
region for H\\b only and without a further classification. Therefore, we have studied 
this low-temperature region in more detail for all three field directions. 

5.1. Measurements in magnetic fields H\\a 

As shown in figure El (a), a positive anomaly of a a (T) at T^ b ~ 7 K signals the Tb 
ordering in zero field. This anomaly shifts to lower temperature as small magnetic 
fields are applied and it is completely suppressed for H > 2 T. The suppression of the 
antiferromagnetic order is in agreement with the ferromagnetic alignment of the Tb 
moments for H > 2 T reported in Ref. [12]. Our data suggest that the Tb ordering 
is of second order since no hysteretic behavior occurs at the phase boundary (see the 
inset of figuraT] (a)). This phase transition can also be detected by magnetostriction 

studies [28]. Magnetic neutron diffraction is still missing. 




Figure 7: Panel (a): FC thermal expansion a a (T) for H\\a showing the suppression of 
the Tb ordering above about 2 T. Panel (b): Magnetostriction measurements Aa(H)/a 
recorded with increasing field below T^ b . 



measurements. The Aa(H)/a curves in figure [7] (b) were obtained with increasing H 
at constant T. Considering Aa(H)/a at 4.5 K, a first step-like expansion occurs at 
H — 0.9 T and a second one at H = 1.7 T. The second one can be attributed to the 
ferromagnetic alignment of the Tb moments [12], but the origin of the first anomaly 
remains unclear thus far. Again, both anomalies show no hysteretic behavior as a 
function of field (not shown). With increasing temperature these anomalies broaden 
and vanish for T > T^ b ~ 7 K. 

5.2. Measurements in magnetic fields H\\b 

For magnetic fields H\\b up to H ~ 1 T, a&(T) shows similar anomalies at T^ h as 
those observed in a a (T); see figures M (a) and [7] (a). In contrast to a a (T), however, 
the antiferromagnetic Tb ordering is not suppressed as H is increased. Instead, another 
ordered phase is induced for H > 1 T. The anomalies in higher fields are sharper 
and significantly larger than those in H < 1 T. For fields H\\b up to 2 T this second 
Tb ordering is stabilized, but with further increasing field the transition temperature 
decreases again and at 4.5 T the transition already occurs below the lower limit of the 
investigated temperature region, see figure [8] (b). The insets of panel (a) and (b) of figure 
M compare measurements of a&(T) obtained with increasing and decreasing T. For both 
types of transitions no hysteresis is observed. In Ref. [12] it has been reported that 
the wave vector of the antiferroamgnetic Tb oredering changes from an incommensurate 
value below 1 T to a commensurate one for higher fields H\\b. Thus we attribute the 
broad low-field anomalies in a&(T) to transitions to the incommensurate phase, while 
the sharper anomalies above 1 T signal transitions to the commensurate phase. 

Panel (c) of figure [8] presents the relative length change Ab(H)/b as a function 
of increasing H. These curves confirm that two different phase transitions have to 
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Figure 8: Panel (a) and (b): FC thermal expansion a&(T) for H\\b showing two different 
types of anomalies at T^ b . Below about 1 T the anomaly is rather broad and hardly 
changes with field. For larger fields the anomaly becomes very sharp and its position 
changes non-monotonically with field. Panel (c): Magnetostriction Ab(H)/b recorded 
with increasing field. For clarity, the curves are offset with respect to each other. 
Panel (d): Magnetostriction Ab(H)/b recorded with increasing and decreasing field 
at T = 4.5 K (see text). 

be distinguished in a&(T) at T^ h . As a function of increasing H, Ab(H)/b shows a 
pronounced step-like contraction at 1.4 T, in agreement with the observed transition 
from an incommensurate to a commensurate wave vector of the Tb ordering [12]. At 
a somewhat larger field we find an anomalous expansion of the b axis suggesting that 
the commensurate phase is left again. With further increasing field another step-like 
contraction takes place, which is due to the LTI-to-LTC transition of the Mn moments 
discussed in section [3T2l The latter is present up to about 25 K (see figure [4]), while 
the transitions at lower fields can be observed only at temperatures up to T ~ 7.25 K. 
Thus, the transition temperature of the commensurate Tb phase is slightly higher than 
TJ b in zero field. 
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T (K) T (K) 

Figure 9: Panel (a): FC thermal expansion a c (T) for H\\c around T^ b . Panel (b) and 
(c): Thermal expansion recorded with increasing (•) and decreasing (-•-) temperature. 

Figure[8](d) shows an expanded view on Ab(H)/b at 4.5 K as a function of increasing 
and decreasing field. As in Aa(H)/a, there is an additional anomaly (of unknown origin) 
marked by ©, before the incommensurate-to-commensurate transition occurs at the 
anomaly marked by ©. The position of the first anomaly is not hysteretic with respect 
to the direction of the field sweep, whereas the position of the second one as well as 
the anomalies around 5 Tesla show some hysteresis. Our phase boundaries qualitatively 
agree with the phase diagram presented by Kimura et al. for H\\b in Ref. [10]. 

5.3. Measurements in magnetic fields H\\c 

For a magnetic field H\\c, we find a complex field dependence of the anomalies of a c (T), 
which signal the various rearrangements of the Tb moments below T^ b . In zero field, 
an anomaly of negative sign indicates the ordering of the Tb sublattice. This anomaly 
broadens as higher fields are applied and can be observed up to ~ 2 T; see figure [9] (a). 
For higher fields, a broad anomaly of positive sign shows up. This anomaly continuously 
increases in magnitude for further increasing H and has been detected up to H = 7 T; 
see also figure [5] (b) . Presumably, the two anomalies of opposite signs belong to two types 
of different Tb phases similar to those observed for H\\b. The need to distinguish two 
anomalies also for H\\c becomes more obvious in panel (b) and (c) of figure[9l There, we 
present the anomalies of a c (T) recorded with increasing and decreasing temperature in 
H = 1.25 T and H = 2 T. For both field strengths, a c (T) shows only one broad anomaly 
of negative sign as a function of increasing T. For decreasing T, however, this anomaly is 
overlapped by a second one of opposite sign. This complex hysteretic behavior prevents 
an exact determination of clearly-defined phase boundaries. Nevertheless, our data 
suggest that the Tb ordering is not suppressed up to 7 T. 



New features in the phase diagram of TbMnOs 



16 



6. Summary 

We have determined the magnetic-field temperature phase diagram of TbMnOs 
by high-resolution thermal expansion and magnetostriction measurements. The 
measurements have been performed in longitudinal magnetic fields applied along all 
three crystallographic axes. The fact that we find rather pronounced anomalies at 
the various phase transitions proves that all these transitions strongly couple to lattice 
degrees of freedom and allows for a detailed investigation of the phase boundaries. Our 
data reveal various new features in the phase diagram of TbMnOs: Firstly, the phase 
boundary between the LTI and the LTC phase for H\\a has been determined for the 
first time. This first-order phase transition shows a broad hysteresis, which is strongly 
enhanced at lower temperatures. Secondly, we find clear evidence that even in high 
magnetic fields H\\a no direct HTI-to-LTC transitions take place, contrary to previous 
reports [10,15]. This is also the case for high fields H\\b. Our data suggest that for both 
field directions the Mn subsystem always transform from the HTI first to the LTI phase 
before the LTC phase is finally established at low temperature. Thirdly, we observe a 
strongly anomalous behavior at the LTC-to-LTI transition for H\\b. Similar to glass 
transitions, double-peak structures show up in ab(T) and 9A q I ^ L when the LTC phase 
is left as a function of temperature and magnetic field, respectively. Nevertheless, we 
could not observe any dependence of these double-peak structures on the sweep rate 
of the temperature or magnetic-field changes during the prior LTI-to-LTC transitions, 
which would be a typical indication for a glass-like transition. 

Besides the positions of the various phase boundaries, our data also yield 
information about their uniaxial pressure dependencies. The uniaxial pressure 
dependencies of the Neel temperature Tjv of TbMnOs have the same signs as those 
of GdMnOs. This confirms our previous conclusions [14], that the increase (decrease) 
of T/v for uniaxial pressure applied along the a (b) axes arises from a decrease (increase) 
of the orthorhombic splitting e, which causes an increase (decrease) of J™- I n order to 
explain the uniaxial pressure dependence for pressure along the c axis, the finite J^ FM 
has to be taken into account. The analysis of the pressure dependencies of Tfe suggests 
that the optimum value of e needed to establish a ferroelectric order with a maximum 
Tfe is located between e of GdMnOs and e of TbMnOs. 

Concerning the ordering of the Tb moments, our data confirm the suppression of 
the incommensurate ordering below T^ b = 7 K for magnetic fields H > 2 T applied 
along the a axis. For H\\b, we also observe that a clear change in the ordering of 
the Tb moments is induced for H > 1 T, in agreement with the incommensurate-to- 
commensurate transition found by neutron scattering [12]. Clear anomalies, due to 
the ordering of the Tb sublattice are also present for H\\c and our data indicate a 
rearrangement of the Tb moments around 3 T. The ordering is not suppressed in fields 
up to 7 T, but a clear attribution to different phases is prevented by the very complex, 
hysteretic field and temperature dependencies. 
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